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We study the evolution of the non-linear curvature perturbation during perturbative reheating, 
and hence how observables evolve to their final values which we may compare against observations. 
Our study includes the evolution of the two trispectrum parameters, ^nl and tnl, as well as the 
scale dependence of both /nl and tnl- In general the evolution is significant and must be taken 
into account, which means that models of multifield inflation cannot be compared to observations 
without specifying how the subsequent reheating takes place. If the trispectrum is large at the end 
of inflation, it normally remains large at the end of reheating. In the classes of models we study, 
it is very hard to generate tnl ^ /nli regardless of the decay rates of the flelds. Similarly, for 
the classes of models in which (;nl — tnl during slow-roll inflation, we flnd this relation typically 
remains valid during reheating. Therefore it is possible to observationally test such classes of 
models without specifying the parameters of reheating, even though the individual observables 
are sensitive to the details of reheating. It is hard to generate an observably large (Jinl however. 
The runnings, Uf^-^ and Wrp^Li tend to satisfy a consistency relation WrNL = (3/2)n/NLi but are in 
general too small to be observed for the class of models considered regardless of reheating timescale. 
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I. INTRODUCTION 



Inflation has become the leading paradigm for solving the horizon, flatness and relic problems in the Standard 
Hot Big Bang picture (for example, see |lH3]) and explaining the origin of structure formation in our universe. 
The original model consists of a scalar field slowly rolling down a flat potential [1,, resulting in an exponential 
expansion of spacetime. More complicated, particle theory motivated models have been studied since then. With a 
vast number of inflationary models in the literature, it is important to constrain and test individual ones in order 
to make connection with particle physics models. Recently, observational constraints on the detailed statistics of 
C, has emerged as a powerful tool for testing different inflationary models. Here C, is the gauge-invariant curvature 
perturbation, quantifying the perturbation in total energy density of the universe. 

Primordial non-gaussianity, as an example, opens up an extra window to constrain different inflationary models. 
While simple single-field models predict negligible levels of non-Gaussianity [H |5] , significant non-gaussianity can 
be generated by different mechanisms, such as features in the inflaton potential [5], the curvaton scenario [Ttillj. 
modulated reheating/preheating P^HTT] , and an inhomogeneous end of inflation [T51[TU]. It is also possible to generate 
significant non-Gaussianity during multi-field inflation 20 23], see [24] . For a complete review of primordial non- 
gaussianity, see |25j . Here we focus on local type non-Gaussianity generated in multifield models via superhorizon 
evolution of the curvature perturbation, C,. 

As emphasised in [2^1 C, continues to evolve after horizon-crossing in the presence of isocurvature modes, and 
therefore so does its statistics. Thus it is important to take into account any superhorizon evolution up to the point 
where all isocurvature modes are exhausted, in order to evaluate the true model predictions for the statistics of C 
to compare against observations. Reheating, as an important part of infiationary model building which involves a 
transfer of energy from the inflaton to the Standard Model particles, may play a role in the evolution of A number 
of previous works in the literature have assumed that reheating is instantaneous such that C, becomes conserved 
immediately |28H31) . This assumption is unrealistic in general as reheating presumably takes finite time to complete. 
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Peterson et.al. [35] have found compact relations between /nl, tnl and (/nl and the tilts of the curvature and 
isocurvature power spectra, under the slow-roll and slow-turn approximation. In particular, they found that for 
detectable non-gaussianity in multifield models without excessive fine-tuning, tnl can not be much larger than f^^. 
These results have also been verified by Elliston et.al. [33] making use of explicit analytic expressions for models with 
separable potentials. Even when reheating is taken into account, we will find that their results continue to hold. 

Using the finite central difference method, we have numerically implemented the SN formalism to follow the 
evolution of C and its statistics for two-field models through a phase of perturbative reheating. In a previous paper |27j , 
we have demonstrated that /nl is in general sensitive to the reheating timescale, whilst is less sensitive, and therefore 
can be considered a better probe of the underlying inflationary model. In general, the sensitivity is model-dependent, 
meaning that single values of /nl can only be reliably used to discriminate between different multifield models if the 
physics of reheating is properly accounted for. Here we extend our previous work to the trispectrum and the running 
of the non-linear parameters n /^^ and fiTNL ■ Besides investigating how reheating changes the inflationary predictions 
of these individual observables, we also explore the possible relations between different observables. The aim is to find 
whether there are consistency relations between observables, which would therefore be a smoking gun of the scenario 
which gives rise to the relation. 

As we found was the case for /nl pTj . we find that the non-linear parameters in the trispectrum are also sensitive 
to the reheating timescale, while ^nl remains too small to be observed in general. The runnings rifj^^ and rir,^^ are 
small except in some cases of the quadratic times exponential model. Though individual observables depend upon 
the reheating timescale, consistency relations between them in some classes of multifield models may survive through 
reheating, offering hope to test such models without specifying the details of reheating. 

The paper is organised as follows: in Section ( [TT| ) and Section ( III ) we introduce some background material including 
the definitions of the non-linear parameters /nlj "^nl and (7nl and the runnings and Uj-j^j^- In Section (IV I we 
recall the 6N formalism and give the formulae of the primordial observables in terms of the SN derivatives. Then, 
in Section (|v]) we study how reheating may alter the infiationary predictions for tnl, 9nl, JT-Znl and Urf^^, including 
relations between them, in canonical two-field models. We study the two classes of models where a minimum exists in 
one or both field directions in Section ( V A I and ( V B ) respectively. Our discussion and conclusions are presented in 
Sections ( VI I and ( VII ) . The reader who is interested in the details of our numerical recipe is referred to our previous 
paper [27] . 



II. BACKGROUND THEORY 



The class of two-field models considered in this paper are described by the inflationary action: 
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where Mp = I/a/SttG is the reduced Planck mass. The standard slow roll parameters are defined as 
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where subscripts denote differentiation with respect to the fields, (pi. The dynamics of the scalar fields are governed 
by the Klein-Gordon equation 



(pi + 3H^i + W^,=0, 



(3) 



where the first term can be neglected during slow-roll inflation. After inflation ends, the fields start to oscillate about 
their minima. If the period of an oscillation is much shorter than the Hubble time, the fields can be interpreted as 
a collection of particles with zero momenta that decay perturbatively to bosons Xb and fermions ipf via interaction 
terms hke — ^g^^p'jxl ^^'^ —hipftpff^- This is the process of perturbative reheating j33] . 

As a phenomenological prescription, we model reheating by adding friction terms Tjipj to the Klein-Gordon equation 
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Eq. (H [27] 

<fi + {iH + Ti)<pi + W^, =0, (4) 

p^ + AHp^ = ^r,^?, (5) 

which couples the scalar fields to an effective radiation fluid p^. The decay terms are 'switched on' only when the 
scalar fields pass through their minima for the first time, with the conditions m^j ^ maxjiif, F} satisfied, where 
r = J2i [SI]- For simplicity, we assume the decay rates Tj are constant and the decay products are relativistic and 
thermalised instantaneously. The completion of reheating is taken to be the time when the universe becomes radiation 
dominated, i.e. fJ-,, ~ 1. At that point, the isocurvatre perturbations have decayed away and become negligible, hence 
observables will freeze into their final values. Very recently a study of the possible survival of an isotropic pressure 
perturbation during reheating has been studied [IS], where the fields are allowed to decay into both radiation and 
matter. Even allowing for this, in all models studied the isocurvature mode quickly became negligible. 



III. THE CURVATURE PERTURBATION, ( 

At leading order, the statistics of C are measured in terms of the power spectrum, bispectrum and trispectrum, 
which are defined in Fourier space by 

(Ck.Ck.) = {27:f5^ik^+k2)^Vc{k^), (6) 
(CkiCk.Ck3) = {27rf6^{ki+k2+ks)B^{ki,k2,h), (7) 

(Cki Ck. Ck3 CkJ = (2^)'<53(ki+k2+k3+k4)T^(fcl,fc2,A:3,fc4,fcl2,fcl3)- (8) 

where kij = |ki + k^ j. Here the delta functions are present due to the assumption of statistical homogeneity and 
isotropy. The level of non-gaussianity can then be parametrized by the dimensionlcss non-linear parameters /nl i ''nl 
and gNL, 

Bc(fci,fc2,fc3) = ^hL[Pdki)Pcik2) + 2perm,s] , (9) 

54 

Tc{ki,k2, fcg, A:4,fci2, /cia) = tnl [P(;(fci2)^c(^i)-Pc(^3) + Uperms] + ^Snl [^c(^i)^c(^2)-Pc('=3) + Sperms] ,(10) 

which are in general functions of wavevectors fc^ and thus are shape dependent. For canonical models however, 
non-gaussianity peaks in the local shape, which is defined by C = Cg + 3/nlCg/5 where Cg is the Gaussian part 
of (. Here we focus on the local shape only, for which the current constraint from WMAP 9-year data on /nl is: 
^3 < /nl < 77 at 95% CL [36]. An even tighter constraint comes from large scale structure —37 < /nl < 25 at 
95% CL [37]. For the trispectrum, WMAP 7-year data gives the following constraints: —0.6 < tnl/10'' < 3.3 and 
-7.4 < ^nl/IO^ < 8.2 [Ml 133, with [ID] finding the compatible constraint -5.4 < 5nl/10^ < 8.6. A slightly tighter 
constraint for ^nl was also found in [37], but with the caveat of how to model the scale dependent bias due to ^nl- 
All of these will be improved considerably by Planck data very soon. In the absence of a detection, the bounds are 
given by |/nl| < 5 gl], tnl < 560 [12], and I^nlI < 1-6 x 10^ [33]. 

Like the power spectrum, it is natural that the non-linearity parameters are scale dependent [3T1 H3H35] , quantified 
by their runnings. For instance, the runnings of /nl and tnl, denoted by Uf^^ and nr^L, are defined by 

_ dln|/NL| 
= ^In^ ' 
_ dln|rNL| 

= -^]^ ' ('2) 

where k marks the length of any one side of the n-gon, provided that all sides are scaled in the same proportion ffl] . 
Examples of models where rif^^ and rir,^^ can be observationally large, i.e. 0(0. 1), are the curvaton model with 
quartic self-interaction terms [Ml [37] and Modulated Reheating [33] . Forecasts have been made to assess our ability 
detect the running of these non-linear parameters. For Uf^^, Planck could reach a 1 — ct sensitivity of f^uf^^ ^0.1 
given /nl = 50 [31]. By measurements of the CMB /x-distortion in a CMB experiment such as PIXIE, /nl and rir^^ 
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could also be measured to an accurancy of the order of 0(0.3) and 0(0.6) respectively for /nl = 20 and tnl = 5000, 
and similarly in large-scale surveys such as Euclid [49j . ^ For any non-linearity parameter, the error bar on its scale 
dependence is approximately inversely proportional to its fiducial value [IS]. 



IV. THE (5N FORMALISM 



The SN formalism |50H52| has been used extensively throughout the literature to compute the primordial curvature 
perturbation and its statistics. The formalism relates C to the difference in the number of e-folds of expansion, SN, 
between different superhorizon patches of the universe, given by |52j (or [531 154j for the covariant approach) 



(13) 



where N is defined as the number of e-folds of expansion from an initial flat hypersurface to a final uniform energy 
density hypersurface. We take the initial time to be Hubble exit during inflation, denoted by t*, and the final time, 
denoted by tc, to be a time deep in the radiation dominated era when reheating has completed. All repeated indices 
are implicitly summed over unless stated otherwise. Here Nj — dN/{d(pi^,), the index / runs over all of the fields, 
and subscript * denotes the values evaluated at horizon-crossing. In general, N{tc,t^,) depends on the fields, fi{t), 
and their time derivatives, (pi(t). However, if the slow-roll conditions, SHipj ~ —Wj, are satisfied at Hubble exit, 
then N depends only on the initial field values. The radiation fluid remains effectively unperturbed at horizon exit 
as it does not yet exist, and so does not feature in the above expansion. 

For canonical models, the non-linear parameters defined in Eq. (9p0 ) are dominated by their shape independent 
parts, which under the 6N formalism are expressed as [SH [S5] 
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Similarly, in terms of 5N derivatives, the runnings Uf-^-^ and t^tnl ^^'^ given by 
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{NmNmT 
-H and the slow-roll field equations. 



NijkWi, = Wjk*-NijWik*-NikWij,-NiWijk* 



(19) 
(20) 
(21) 



where Eqs. ( 20|2T ) are derived by differentiating Eq. (19) with respect to Lp*j. The results that during slow-roll, higher 
order 5N derivatives can be eliminated in favour of lower order ones whenever they come in combinations with ^pl^, 
such as Nitfi^,, Njjifii^, was first noticed by Lyth and Riotto, for instance see Eqs. (113) and (114) in [57] where they 
have used these to work out alternative expressions for and its running. Following a similar approach here, we 
n /pjL and n7-j,L ■ This allows us to rewrite n f-^-^ and ni-^L in terms of only first and second order 



extend it to the case of 



^ Their definition of n^p^L differs from the one used here, in the fact that in their case, rirj^^ ^ din |TNL|/(<ilii fc). The two definitions are 
related when the four k vectors form a square by 2?!'^''^^'''' = n'^^, in which case we have to double their forecasted error bars when 
comparing to our definition of rir^L • 
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derivatives of N as follows 
'^/nl= -2K-l + 2e,) 

riTNL = -3("C - 1 + 2e*) 
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(23) 



Here n,^ — 1 is the spectral index. Eqs. ( 22p3 ) are two useful results of this paper. Whilst Eqs. (22 23) are equivalent 
to Eqs. ( [iTpSl ), they possess significant computational advantages over the former since they involve lower order SN 
derivatives which are relatively easier to evaluate in general compared to higher order ones. 



V. SENSITIVITY TO REHEATING 



In this section we present numerical results for the statistics of ^ for the class of two-field models where a minimum 
exists in one or both field directions, focusing on those models which can produce large values of /nl and tnl during 
inflation. We focus on the trispectrum and the running of /nl and tnl. In what follows, x may be identified as the 
inflaton and ip as the field which sources the isocurvature perturbations. For the one minimum case, the ip field is 
not directly involved in the reheating phase and so F;^ = at all times. For models with two minima, both fields can 
decay to radiation and so both F^ and T^p can be non-zero. 



A. One minimum 



First we consider a two-field model where a minimum exists in only one of the field directions. In particular, we 
study the 'runaway' type model 

W{v,x) = Wox'e-^'^"/K . (24) 

This model was first introduced in [55] and has been studied extensively in the literature since then [551 [55M5] . It does 
not possess a focussing region where the bundle of inflationary trajectories may converge, meaning the isocurvature 
mode would never be exhausted unless reheating is taken into account. By placing the ip field close to the top of 
the ridge at horizon-crossing, a large negative /nl can be produced [55] . Here we restrict ourselves to the parameter 
space where a detectable level of non-gaussianity, |/nl| > 0{1), is generated by the end of infiation. ^ 



1. The trispectrum 



Before studying how the non-linear parameters tnl and (/nl evolve during reheating, it is useful to consider their 
evolution during the inflationary phase. Because the potential is of product-separable form, analytic expressions exist 
for Tnl and ^nl during slow-roll. These have been studied extensively in [33]. Anderson et. al. [60] have also studied 
the evolution of tnl and (7nl in this model using the moment transport equations developed in [63) . To summarise, 
a large tnl is produced in similar regions of parameter space as that of a large /nl- 5nl remains subdominant 
throughout inflation except possibly if there are significant terms beyond quadratic order in the potential. 

Here we are interested in the post-inflationary evolution during reheating. In particular we study how tnl and (7nl 
evolve with different decay rates F^^,, and how their final values at the end of reheating depend on F^^. We start with 
Tnl- In Fig. ^ we show the evolution of tnl during reheating for two different decay rates F^, for two slightly different 
slopes of the ridge in the potential as determined by A. The model parameters arc A = {0.05, 0.06}, (p^, — 10"'^ Mp 



^ We have also studied two slightly different models where significant terms beyond quadratic order are present. The potentials are given 
hy W = Wox'^e~^f /*^p and W = Wox^e~^^ /*^p. The qualitative behaviour for these models is similar, with jnl negligible in 
general. 
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Left panel: The evolution of tnl during the post-inflationary period, with 
IG.OMp. Right panel: Same initial conditions with A = 0.06. The solid vertical line denotes 



FIG. I: Potential: W{<p,x) = WoX 
X = 0.05, ip* = lO-^Mp and ; 

the end of inflation, A'e, and the dashed line denotes the start of reheating, N\^=o here and in all subsequent figures for this 



one minimum model, 
infiatiori value. 



Note the final value of tnl can either grow or decay with larger T^, and are different from the end of 



and X — 16Afp. As we found with /nl in |'^7j . tnl oscillates during reheating when x oscillates about its minimum. 
No generic trend independent of A can be seen as the decay rate is increased, as we see that the final value of tnl 
can either grow or decay as the reheating timescale increases. This may be understood by making approximations in 



Eq. ( 15 ) and determining how the SN derivatives evolve as follows: 

As we demonstrated in [2^, during reheating the N^^ and N^^ are negligible compared to iV^^. Together with the 
scaling relation found between N^^ and N,p, where N^^ w N^/ip^, |27], we may then write tnl as 

Here = = (2e* )~^/^, which for this potential, is independent of A. The result that w gi^ = const comes from 
the fact that the x field dominates the energy density over the whole evolution. This algebraic function has three 
stationary points at certain values of N^, 

N^^O,±y^. (26) 



The Nip ~ root is an inflection point where tnl = 0, while the iVi^ = ±Y^2g^ corresponds to a local maximum. Both 
— and N^p = \f2gl roots are unphysical here because A^^ is always negative with diverging trajectories. The 
other root, A^;^ = — y^2g^, however is physical and bounds the maximum value of tnl, given by 



This bound depends entirely on the initial conditions at horizon-crossing, not on any superhorizon evolution including 
reheating. The final value of tnl at the end of reheating is of course dependent on however. But since a bound 
exists, even if the details of reheating such as F;^ are unknown, we are still able to constrain the possible range where 
Tnl could lie in this model. 

We now repeat the same analysis for (7nl- In Fig. [2] we show the evolution of ^nl during reheating for two different 
F;^. The model parameters are A = 0.06, (/s* = 10 Mp and x — 16Mp. Similar to /nl and tnl, we find that ^nl 
oscillates during reheating, with the final value at the end of reheating sensitive to the decay rate T^. Just as for the 
case of the second order 5N derivatives, we have found that there is also a hierachy for the third order 5N derivatives 
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FIG. 2: Potential: W{ip,x) = Wox^e ^"^ /*^p . The post-inflationary evolution of jnl for three different decay rate F^. The 
model parameters are A = 0.06, (p, = 10"^Mp and x* = 16.0Mp. Left Panel: A = 0.05, Right Panel: A = 0.06. Compared to 
tnl, 5nl remains subdominate after reheating and is consistent with being zero in future experiments, though it differs slightly 
from the end of inflation value as shown in Table IH 



with N.r, 



being much larger than the other third derivatives. Using this, Eq. ( 16 1 can be reduced to 



5nl 



25 N^vf^v 
54 {m + 7V2)3 



(28) 



Compared to tnl, however, it remains subdominate (<C O(IO^)) and too small be observed in ongoing CMB experi- 
ments. 

2. Runnings of non-linear parameters, n fj^-^ and Ur^j^ 

Next we study the runnings of the non-linear parameters, rifj^-^ and Tirj^j^- In Fig. [3] we give the whole evolution of 
n/jjL 9.nd n-Tj^j^ including reheating. The model parameters are A — {0.05, 0.06}, ip^, = 10~^Afp and x — 16Afp, for 
the decay rate = \/10~^WoMp. For A — 0.06, we find that both n/„L ^-^^ rir-f^-^ are too small to be observationally 
relevant for CMB experiments, regardless of the decay rates F^. For A = 0.05, however, n/^j^ and r^TNL are much 
larger and of order 0(0. 1), which could be potentially observed by Planck provided that the fiducial values of the 
non-linearity parameters are large enough. 



To understand why n/^^ and n^-^L are much larger for A — 0.05, we first rewr ite Eqs. ([22][23| as 
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(29) 



, (30) 



where r is the tensor-to-scalar ratio 



From this, it is not difficult to see that the second terms 



usmg r = ^^^^^ 

in the first fine of both equations are small in general if r satisfies observational constraints. Making use of the 
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FIG. 3: Potential: W{ip,x) = Wox'e-^'^ '^'^ . Left panel: The evolution of JT-z^l- Right panel: The evolution of Ur^-^. The 
model parameters are A = {0.05, 0.06}, ip, = 10~^Mp and x* = 16.0Afp, and — VlO'^WoMp. For A — 0.05, n/j^^ ^-nd Ut^-^ 
may be large enough to be observationally relevant, while for A = 0.06 the non-linear parameters are almost scale-independent. 




FIG. 4; Potential: W{ip,x) = Wox^e"^*^ /"p. Comparison of the exact Eqs. 
panel: The evolution of n/^^^. Right panel: The evolution of Wtnl- The model parameters are A — 0.05, (p„ = 10~^Mp and 
X* = 16.0Mp, for the decay rate = y/TO~^WoMp. The two equations agree to a good approximation after about 30 e-folds 
of inflation. 



29I3O1 and approximate formula Eq. |3l|. Left 



approximate formulae for /nl and tnl, the dominant terms in Eqs. ( 29]|30 ) are 

3 



X 



(31) 



We have also assumed that the numerators in the square brackets in Eqs. ( 29p0| are dominated by the and 
terms. In Fig. |4j we show the comparison between the exact Eqs. ([29|30) ana the approximate formula Eq. (31 1. 
From this, we can see the approximate formula agrees very well with the full expressions after about 30 e-folds of 
inflation, even during the reheating phase. From Eq. (31), one may see that the runnings are relatively large when 

which is the case when A = 0.06. 



N^, which is the case when A 



0.05, but very smaU when |iV<^| > |A^;^ 



Notice that if \N^ \ 3> \N^\, the runnings are driven to zero and hence become independent of the decay rate. 
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FIG. 5: Potential: W{x,v) 

model parameters are A — {0.05, 0.06}, ip* — 10"'^ Mp, x 
about 30 e-folds of inflation, showing the consistency relation Eq. ( 32 1 is satisfied. 



p. The evolution of the ratio nr^^/nf-^-^ until the completion of reheating. The 

= V0.3WoMp. The ratio settles to 3/2 quickly after 



16Mp and F;^ 



Whether n/^j^ and tIt^-^^ are of a detectable level or not, we find that they satisfy the consistency relation 

3 

regardless of F^^ and thus the reheating timescale. This relation was first observed to hold for some classes of two field 
models in |44) . We provide an example of this scaling behaviour in Fig. |5] where we observe it to hold throughout 
most of the evolution, except partly during the first 20 e-folds (of course this evolution is not itself observable, only 
the final values). We will discuss this relation in further details in Section [Vl| 

In Table |l] we summarise the results, showing the comparison between the primordial observables evaluated 
at the end of inflation and at the end of reheating. This is one of the main results of this paper. Notice the different 
qualitative behaviour for the non-linear parameters in the models for different A, where the magnitudes of /nl and 
tnl decrease with larger F^ for A = 0.05, but increase for A = 0.06. In general, the final values of the non-linear 
parameters at the completion of reheating is different from the end of inflation values, whilst g^i^ remains small 
<SC O(IOO) in this model which is unlikely be observable in future experiments. The runnings Uf^-^ and j^tnl are large 
in the case A = 0.05 and are redder for larger F^. 
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ing,A = 


0.05 






End of Reheating, A = O.Of 








/nl 


T'NL 


5nl 








/nl 


■''NL 


5nl 




"tnl 




Vio-5 


-33.4 


2.25 X 10^ 


-13 


-0.105 


-0.157 


Vio-5 


-4.35 


28.1 


-2.41 


-9.1 X 10-'' 


-1.3 X 10- 


-3 


Vio-3 


-31.5 


2.27 X 10^ 


-11.6 


-0.137 


-0.205 


Vio-3 


-5.54 


44.5 


-2.62 


-1.4 X 10^3 


-2.1 X 10- 


-3 


Vio-1 


-26.9 


2.01 X 10^ 


-9.96 


-0.177 


-0.266 


Vio-1 


-7.14 


73.9 


-2.96 


-2.3 X 10-3 


-3.4 X 10- 


-3 



TABLE I: Statistics of C for W{lp,x) = Wax^e'^'^ ^'^'p for difi'erent decay rates. AU decay rates are in units of -/WoM^. We 
give values computed at the end of inflation {t^) and at the completion of reheating (final) where is conserved. The model 
parameters are A — 0.05 (Left panel) and 0.06 (Right panel), if* = 10"'^ Mp and x* ~ 16.0A^p. 
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FIG. 6: Potential: W{(p,x) ~ Wo[\rn?x^ + A*(l — cos(y^(/3))]. The evolution of tnl during post-inflationary period. The 
model parameters are A* = m^/^/47r^, f — m — 1, ift, = — 0.001)Mp and x* ~ 16.0Mp. For these model parameters, the x 
field minimises before the ip field. Left panel: Equal decay rates, — T^p. Right panel: Unequal decay rates, F^^ ^ F^. The 
solid vertical line denotes the end of inflation, A'^e, and the dashed lines denote the start of reheating, N\^^o (blue) and A'^|x=o 
(black), respectively, in this figure and all subsequent figures for this two minima model. Notice that tnl changes by a few 
orders of magnitude during reheating. Also, tnl is sensitive to F^ and F^ if there is a hierarchy between the two decay rates. 

B. Two Minima 

Next we consider a model where the potential has minima in both field directions. Both fields can decay to the 
eff'ective radiation fluid and are directly involved in reheating. The potential considered is the efi^ective two-fleld 
description of axion N-flation, where the potential given by |26j 



W{ip,x) = Wo 



m^X^ + A"' 1 - cos 



27r 

7 



(33) 



The axion ip, is described by its decay constant / and its potential energy scale A**. To generate a large non- 
gaussianity, we must have ip close to the "hilltop" at horizon-crossing 26J. In this configuration, the second field x, 
drives infiation."^ 



1. Trispectrum 



The model parameters we consider are A4 



^2/V47^^ = (I - 0.001)Mp, x* = 16Mp and / = m = 1. Each of 
/nl, Tnl and (7nl are negligible during inflation as the axion if is sufficiently light that it remains almost frozen near 
the top of the ridge. In this sense, this scenario is similar to the curvaton model. Things are different after inflation 
ends however. When inflation ends, the axion ip starts rolling down the ridge, producing a negative spike in /nl- 
/nl then evolves to positive value when the p field converges to its minimum |27l . It is similar for tnl, except tnl is 
always positive. In Fig. [6] we give the evolution of tnl during reheating for various combinations of and T^p. 

Similar to /nl as studied in 27J, we find that although the final value of tnl is different from that at the end of 
inflation, it is almost completely insensitive to the decay rates if T^ — T^. Things are different however if there is a 
mild hierachy between and r,^. When T^ 7^ Tip, the final value of tnl does depend on the reheating timescale. 
Compared to the value where = T^, it grows for F^ > F^ and decays for F^^ > F^. 



^ We have also studied the models W - 
conclusions were found in these models. 



Wo [|m2x4 + A4 (l-cos(2l^))] 



and W - 



Wo (jr9X* + Vo + hip+ |f A(/33 + ^jiip*). Similar 
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Unlike the previous model in Section V A no scaling relation is found between N^^ and N^. Yet we can still make 



use of the observations that N^^ and N^, dominate over Nyy,Nu,^ and Ny respectively to rewrite Eq. (151 as 



Ar2 



(34) 



For 5nl, things are similar to /nl and tnl- In Fig.jTjwe give the evolution of ^nl for different combinations of and 
r^, with the same model parameters. While the final values of ^nl at the end of reheating are different from that at 
the end of inflation, they are almost completely insensitive to and F^ unless there is a mild hierachy between the 
decay rates. Again with a hierachy between the third order 5N derivatives found, ^nl can be well approximated by 




N 



FIG. 7: Potential: Wix^'P) = H^o | |m^x^ + A** [l - cos(^<y9)J |. The model parameters are A* = m^p/An^, 99. = (| - 
0.001)Afp, X* = 16Mp, f = m = 1. Left panel: Equal decay rates, F^ = F^; Right panel: Unequal decay rates, F;^ 7^ F^. 
Similar to tnl, <?nl changes by a few orders of magnitude during reheating and is more sensitive to the decay rates whenever 
there is a hierarchy between them. 



25 -/Vyyy . 

^^"^54^- ^^^^ 

During slow-roll, EUiston et.al. |33j have shown that ^nl is roughly of the same order as tnl and the following relation 
holds 

^5jvl ~ Tnl , (36) 

for non-vacuum dominated sum-separable potentials, given that tnl is large. Here we have found that this holds 
beyond the slow-roll regime and during reheating for a wider range of mass ratios between the axion and inflaton 
where they both minimise after the end of inflation. For example, from FigjSj we can see that this relationship is only 
mildly violated when F^^, ^ F,^. 



2. Runnings of non-linear parameters, n and Wrp^L 

We now turn our attention to the study of the runnings nf^-^ and r^rNL i'^ ^^i^ model. Similar results are found 
as in the one minimum case where A = 0.06, i.e. both runnings are small and — '^tnli except the relation 



Eq. (32) may be mildly violated when there is a hierachy between the decay rates. 



In Table [IT] we summarise these results, showing the comparison between the primordial observables evaluated at 
the end of inflation and at the end of reheating. This is one of the main results of this paper, which clearly show the 
non-linear parameters in multifield models strongly depend on the reheating timescale in general. Notice the large 
differences between the statistics evaluated at the end of inflation, compared to the end of reheating. This is because 
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r =Vio"^, r = Vio"* 
r =Vio"'^, r =Vio"^ 
r =r =Vio"'^ 




67 
N 



,r = ^10-" 



r =^10"*, r =V10"' 

- y. 

r =r =^10"" 



69 



70 



FIG. 8: Potential: W{x,v) = Wollm^x^ + A** [l - 008(^1/?)] }■ The model parameters are A" = m^f/Ait^, 93. = (| - 
0.001)Mp, X* ~ 16Mp, f = m = 1. Left panel: The evolution of the ratio (54/50) ((;nl/tnl) during reheating for different 
combinations of decay rates. Right panel: The evolution of the ratio ?^tnl/'^/nl during the post-inflationary period. Notice 
that the relations Eqs. ( 32 1 and ( 36 1 are satisfied after reheating, only mildly violated when ^ 



the axion field only begins to roll after inflation has ended and so until this point, the observables do not evolve 
appreciably. 



End of Inflation 







./nl 




5nl 














5.9 X 10^3 


1.3 X 10-3 


-3.1 X 10-5 


1.7 X 10- 


2 


4.9 X 10" 


4 












End of Reheating 














/nl 


•''NL 


5nl 
















6.88 


0.69 X 10^ 


0.63 X 102 


-1.2 X 10 


-6 


-1.8 X 10 


-6 


Vio- 


^2 


Vio- 


-2 


6.59 


0.76 X 10^ 


0.55 X 10^ 


-9.3 X 10 


-7 


-1.6 X 10 


-6 


Vio- 


-2 


Vio- 


-4 


4.37 


0.29 X 10^ 


0.29 X 10^ 


-7.2 X 10 


-7 


-1.2 X 10 


-6 


Vio- 


-4 


Vio- 


^2 


13.66 


2.75 X 10^ 


1.91 X 10^ 


-2.5 X 10 


-6 


-3.7 X 10 


-6 


Vio- 


-4 


Vio- 


-4 


6.83 


0.68 X 10^ 


0.59 X 10^ 


-1.1 X 10 


-6 


-1.7 X 10 


-6 



TABLE II: Statistics of ( for W{ip, x) = |^|"i^X^ + A** ^1 — cos (7^'''))] ^"'^ different decay rates. AU decay rates are in 

units of "v/WjA/p. We give values computed at the end of inflation (te) and at the completion of reheating (final) where C, is 
conserved. The model parameters are A* = f'^ / Atv'^ , ipt. — — 0.001)Afp, x* = 16Mp, / = m = 1. Note that the values in 
the second row where F;^ = F^ = do not correspond to end of reheating since the decay rates are zero. However an adiabatic 
limit is still reached as both ip and x behave as matter fluids when oscillating about their minima. 



VI. DISCUSSION 



A. Relation between tnl and /nl 

In general, TNL(fci, ^2, ^3, ^4) and /nl(^i, ^2, ^3) are functions of external momenta which cannot be compared 
directly. Yet in canonical models, when the non-gaussianity is large, it is dominated by the shape independent parts. 
It is thus reasonable to compare the non-linear parameters directly in such models. 

The Suyama-Yamaguchi inequality [16j . for instance, relates /nl in the squeezed limit (fci — > 0) to tnl in the 
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collapsed limit (fci + k2 ^ 0) 

TNL > (^/nl)' . (37) 

This ine quality has been studied and verified extensively in the literature, see e.g. |65H71j . While the equality in 
Eq. (37) holds for single-source models [B5j, multifield models in general give tnl > (6/nl/5)^. For a recent review, 



see [221 

Recently, Peterson et.al.[35] have shown that tnl is not much larger than /^^ i^i two-field canonical models in 
general by applying both the slow-roll and slow-turn approximations, except in cases of excessive fine-tuning. It was 
later shown by EUiston et.al. [33] to hold also for separable potentials. Here we find that this remains true even if the 
subsequent reheating phase is taken into account. 

For the canonical models studied, a large non-gaussianity is typically generated by having a field rolling down an 
extreme point like a ridge or a valley. In these cases, /nl and tnl in general can be approximated by 



6 N^^K 



2 



5 + ^2)2 

{m + iV2)3 



From Eqs. ( [38p9l ), we then have 



1 + ^- (40) 



(6/5)V5l 

As a result, in order to have tnl ^ /nl I/nl| > 0(1), one typcially needs \N^p \ <^ \N^ \ while \N^^\ ^ l^xxl' l-^vxl- 
These conditions are highly non-trivial for any function of N, and in general is difficult to accommodate in canonical 
two-field models. 



B. The relation between tnl and gNL 

For non-vacuum dominated sum-separable potentials, by making use of the analytic expressions for SN derivatives, 
EUiston et.al. ^33, have shown that 5nl and tnl are of the same order 

27 

^9nl » Tnl , (41) 

in the absence of significant terms beyond quadratic order in the potential. The effective N-flation model we studied 
in Section |VB| is of this type. We found that this relation holds not only during the slow-roll regime but also after 
reheating, given tnl is large. The reason that ^nl ^ tnl regardless of subsequent evolution beyond slow-roll may 
be understood if we split the contributions to the non-linear parameters into instrinsic terms which depend on the 
instrinsic non-gaussianity of Sipj at late times and gauge terms which do not. This could be seen in the moment 
transport techniques developed by Mulryne et.al. |63] . where ^ is evaluated by evolving the field correlation functions 
from horizon-crossing to the time of interest, then gauge-transforming to C on uniform energy hypersurface, as in (73) . 
In particular, one can see that the second terms in the moment transport expressions Eqs. (61) and (62) in |60j for 
Tnl and ^nl would be of the same form up to some numerical factor of order 0(1) if they are dominated by one of 
the field bispectrum contributions, i.e. (SipdipSip). Thus it would be expected that the consistency relation ^nl ^ tnl 
should hold so long as the second terms both dominate in the full expressions for tnl and 5nl, even though tnl and 
(7nl may still evolve in time. We intend to return to this in the future. 



C. Comments on ijnl 

So far for all two-field models considered in the lierature, ^nl is at most of the same order of magnitude as tnl and 
is much less than the current observational limit in CMB experiments and large scale surveys which is about O(IO^). 
It has been already shown recently by EUiston et.al. [33] in the analytic slow-roll expressions for ^nl in separable 
models that it is hard to engineer a model where i^nl can be as large as O(IO^) during inflation and dominates the 
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statistics in the trispectrum, even if one goes beyond quadratic order in the potential. The reasons for this can be 
summarised as follows: 

Looking at the analytic expression for (^nLj given in Eqs (3.10) and (3.12) of 33J, all the terms are multiplied by 
second order slow-roll parameters which are of order O(10~'*) in general. In order to have large (Jnl, we need the 
prefactors and gi to be much larger than 10, in particular > O(IO^) if we want ^nl ~ O(IO^). Yet some of the 
prefactors and gi are bounded from above by 10, and even for those which are not, it requires extreme excessive 
fine-tuning for them to be of order > O(IO^) compared to the conditions required for having large observable /nl £^nd 
tnl- Moreover, the region of parameter space where that extreme fine-tuned conditions can be realized in general 
coincides with those where quantum fluctuations become important over the classical drift of potential flow. As a 
result, it is difficult to engineer a model where (7nl ~ O(IO^) during slow-roll inflation in multifield models. 

It remains to be seen beyond the slow-roll regime though. In particular, (/nl could be dramatically enhanced 
such that it is above the observational limit after reheating. Yet we have found the same conclusion here even with 
reheating taken into account for all the models we studied. In some cases, ^nl does increase dramatically from to 
O(IOO) for some of the combinations of decay rates, for instance see Fig. [6] for the axion model. It may be expected 
that a larger hierarchy between the decay rates may thus produce a large observable (?nl- However this is beyond the 
current numerical capabilities of our code. 



D. Relation between rif^^ and n 



The consistency relation Eq. ( 32 1 follows from the class of two-field local type models with C of the form 



C(fc) = C^-^ + Ck^"" + UiC""'^ * C^^^k + gUC"'^ * C°^^ * C^^^k , (42) 

when and g^p are scale independent and C*^'*^: C'^''^ ^"^^ Gaussian variables. For all cases we study only one of the 
fields has a significant non-Gaussianity, so they may fit this ansatz. The question is whether and g^p are scale 
independent, for the models we study. They are if the field which generates non-Gaussianity is strongly subdominant, 
has no interactions with the infiaton field and a quadratic potential. Many of the models we study are approximately 
of this type, and hence we often observe 3n — 2fiTNL • 

For single source models there is a different consistency relation, which trivially follows from tnl oc /^l, 

"tnl = 2n/„L • (43) 
In the limit that (9'''^ <^ Ck'^^ which corresponds to iV^ ^ N'^ ^ 1, the model becomes effectively single source. If 



the assumptions related to (42 1 remain valid, the non-linearity parameters have to be scale independent. 



VII. CONCLUSIONS 



We have studied the evolution of the curvature perturbation through reheating, for the first time going up to third 
order in perturbation theory. This allows us to study the evolution of several observables during this period for 
the first time, namely the trispectrum consisting of two non-linearity parameters, ^nl and tnl as well as the scale 
dependence of /nl and tnl- The calculation during reheating is complex and requires numerical techniques, which 
to date has led to this field being rather neglected. However it is clearly very important, since reheating is required 
after inflation and observables will often not have reached their final value during inflation. It is of course only the 
final value which we may compare to observations, the evolution before this time is unobservable. In our reheating 
model, in which all scalar fields decay into radiation, the isocurvature mode will necessarily decay during this time 
and the curvature perturbation is thereafter conserved. 

Of course the isocurvature mode does not have to decay during reheating in all models, for example it could be 
sustained by giving the infiaton fields multiple decay channels. But allowing them to decay into both matter and 
radiation alone does not appear to stop them decaying quickly [351 . Even if there was further evolution after reheating, 
our work is not redundant, however one would have to continue calculating the evolution until a later time ^. 



* Note that in addition to the uncertainty of how observables are influenced by reheating, there is also an intrinsic uncertainty between 
the predicted global values of observables, and those which we measure in our Hubble volume. This affect is especially strong in models 
with local non-Gaussianity, due to the coupling between long and short wavelength modes |74H76I . 
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As we have found in [57] previously /nl (see also [21]), we find the trispectrum will in general be sensitive to the 
decay rates during reheating, although in some cases in which both fields oscillate during inflation, the sensitivity 
to the decay rates can be very small provided that they are equal for both fields. In general the evolution during 
reheating is large enough that a comparison between observables and their values at the end of inflation would lead to 
the wrong conclusions, since the change in observables may be larger than the expected error bars of the observables. 
While the evolution to an adiabatic attractor during inflation often (but by no means always) results in a small 
value of /nl [SSHZZHZH, this is not the case during reheating, typically a model which is non-Gaussian at the end of 
inflation will remain non-Gaussian, and in most cases which we studied, the sign of the non-linearity parameters will 
also remain the same. The reverse is not always true, we have seen how in the axion model the perturbations are 
Gaussian at the end of inflation but not at the end of reheating. It would be interesting to study how generic these 
conclusions about the survival of non-Gaussianity are. It would also be interesting to study more realistic models of 
reheating, to include a period of non-perturbative preheating. However such studies are very difficult even following 
single field inflation, and typically require a lattice simulation, which goes beyond the scope of this work. 

Despite the evolution of all observables during reheating, we may still hope to test models of inflation against the 
new observational data. We have previously shown that typically the spectral index, — 1 is a more robust observable 
than non-Gaussianity, since it tends to be less sensitive to the details of reheating. With the non-Gaussianity, we 
may instead look for consistency relations between the five observables which we have studied. First of all we have 
studied the well known Suyama-Yamaguchi inequality, tnl ^ (6/nl/5)^, more specifically how strongly the equality 
may be broken. We have found that it is very hard to generate tnl ^ /nl' consistent with other studies. We have 
provided some analytical and quite general insight into why this is. Interestingly, we have also found that for the two 
field inflation cases in which [33j found (/nl — ''"nl, that this relation typically remains true during reheating. Given 
the observational bounds on tnl, it will be hard to observe ^nl in such models. In fact we have not found any models 
with very large (7nl, despite studying several examples of models with a strong self interaction, so that their potential 
is far from quadratic. Finally we have also observed the relation between n/„L arid Ht-nl, showing that in many cases 
3ri/„L — 2n.rpjL both during and after inflation. These relations between observables allow the model behind it to be 
tested even when one cannot predict the actual values of any of the individual parameters. 
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